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ABSfBACT 

ttie  object  of  this  thesis  is  thi-et  fold,    to  present 
the  theoretical  i-iii^:iogy    ipttt-ct-n  the  fluw   oi'   a  coEii)rei>sxDj.fe 
g&s,   such  as  air,   £.iid   the  flow  of  writer  with  &  free  surf&ce, 
to  describe  tht;  details  of  coRstruction  of  th©  wiXbt  t^ble 
which  utilises  the  above  analogy,  and  to  otitlliie  aial  deaoxi- 
strate  the  applicatlcm^  of  the  water  table  in  the  iiiTestl- 
gat ion  of  air  flow* 

The  analogy  betweeo  air  flow  aafi  sf&ter  isith  &  fre« 
surface  is  presented   iia  Part  I.   Fart  II  con  tains  a  dc^scrip- 
tion  in  detail  of  the  cotistruction  of  the  water  table  at 
the  Rei3;sael«er  IPolytechaaic  Institute*   The  applications  of 
the  »eter  table  and  one  dcasKJnstration  of  its  use  are  found 
in  Part  III. 

This  work  ?.as  done  during  the  months  of  febriaarji 
through  !^ay  1947  at  the  E^assel&er  Poiytechnie  Institute. 


th«  follc^wlBf  iKc»a«mcIc  tur«  Is  used  tn  this  p^p^rt 


Cv  »  a^#elfie  h<?'.-.i.  .tat  can:» '^i-tiK..  voXusw^  Btu  i,"er  lb  p#r  di^  F 

d  -  totiil  differential 

d  *-  partial  fiifr^^rt^tlftl 

g  «j  .iiecel@rativ.=i:  .ji   ;^.t:.  vit^r j   <:.  .  *^,    i*p«  p«r  se# 

J  sr  7?8»a5  ft-Xb  ^^r  Btti 

It  -  n^chsmlci^l  work^  Btu  per  Ife 

Q  rr  heat  i^da^id,  Btu  per  lb 

E  s  ga«  coTiStaati    5S«3  for  4ir 

7  r  absolut€  static  t^Bp^r&ture,  d&g  ]^«eikl&tt 

¥  ST  specific  V'oiiaMi  ee  ft  per  lb 

■?  =  velocity,   f^^ 

▼g  =  velocity  of  a  scyusdi  wsire,  f|i# 

^»a3t  r  a^^^ua  velocity 

W  K  weight  floii,  lb  per  s&e 

y  »  ratio  of  $i^  »p«elfie  h^^ts,  C^/Cj^  I  1.53£?S  for  «oM  air 

^  »  SsasK  dmisity,  »lugs   pur  &u  ft 


MmmctAtmB 

*»^#«  =  ccmpcmeiftts  of  the  velocity  Ui  the  ^^T*^  fiirectlor*« 

h  *  water  de»>th,  ft 

Uq  z  tov-.. i  ii<c-..xi^   *,*.r;. ver  depth  -h^irj:  v"     ^,    .t 

feo  r  totiii  h€?::>4  after  igrar^-tiiXie  J^ai^f  ft 

F  r  pressure^  pisfa 

Q  -  qtxsntity  Df  flo?;,  cti  ft  pBw  s%^ 

▼  s  velocity,  fps 

▼*  s  velocity  of  pro|»ag£tio«i  of  a  surf^e®  %'*tve,  fps 

▼*.«^  r  »&ximu»  velocity 

^p7t^  r  re^ta&gular  eoordln&tcs  i»  flow  sp&e« 

^F^T  r  stagnation  8t«^t«  (^ii^seript  ^o*} 

(^  m  sitss  d€nsitgr#  Bli2gs  |Nir  en  ft 


i-t 


fhe  opera tiari  of  &   super i- or;!:;  Viind  tunnel  for  test 
work  at  trsjassonio  and  stipersonte  air  sp©«d«  Is  an  t^xpenslv© 
project,  re  Hiring  costly  e  "saipsent,   elEborste  iristruments- 
tiorif  axid  en  enorisuus  asoutii.  oi  i.   a  i;'i^eabj.a  3zi.ii  of 

skilled  techRiclaBS  and  much  tliB»  ju"«  ne«d@d  isi  the  ft^brlca- 
tlon  of  aodels  to  be  used  in  the  timrif:!  ?nd  In  the  actual 
operation  aiKi  B&intaiance  of  the  txxrai'^j.  .  isa  its  adjimet** 

?fe«  reduction  in  the  amount  of  test  wcxrk  that  sm^t  b* 
done  in  high  veloeity  air  would  save  i^   conslGerahl©  e^iotjit 
of  tiisie  and  expense.  This  would  facilitatfc  design,  experi- 
mentation, and  theoretic«^l  resf^arch  in  th#  field  of  high 
Telocity  airflow. 

It  is  the  intent  of  this  paper  to  &hoff  th»^t  the  rater 
table  say  be  sdvant-.geously  eaplofyed  In  eonjunetion  ?^ith 
t^e  wind  tunnel*   Its  use  for  prelininsry  nork  vould  repl&ce, 
to  a  considerable  extent,   thf^    >  ^-^-   *long  aecessaj'y  in  the 
wind  turmiil. 

The  feork  herein  presented  was  tsndertAken  with  ttiree 
lirlsary  objectives  in  mirdi    (1)   ^r.  investig.ntion  of  the 
existing  theoretical  analogy  bet  ten  the  two  ciiiaen^xon^x 
flow  of  a  e(»[ipressible  gas,  sueh  as  sir,  and  the  flois  of 
«hallo-r  ^rter  Kith  n     free  sixrf.?ee.    f^>   the  constru©ti<m 
of  a     T-.?Ater  table ,   ^  'at^YXc-  wo  utj.xxj&e  the  flow 

aaalogy,   and   (5)   an  i  tion  of  the  flow  around  ft 

few  b^.sie  isodels  in  or^er  to  check  the  smKlogy  with 
existing  air  iiov    x,hfeory  and   to  de»on»tr&te  the  uscfulnes* 
of  the  water  table. 


W  ^  coaparisoxi  of  the  energy  atn^  coritinuity  equations 
tor  two  ^Limmoator'^l  gas  flo^^?  $.n6  water  flar  rith  z.  free  sur- 
face,  it  a-as  fouiiu    uiiw,i   '.ii  -jaialogy  ^^^^sus   betsefiaa  tn©  r&Lio 
of  absoiutis  air  tejEBr?erature»  sjad  ii«tt^r  deptk  ratio,   th© 
r&tio  of  gas  des^ltieft  aiid  W4ter  depth  ratio,  ^icd  thB  ratio 
of  ga.s  p-ressiares  aad  th^  squsre  of  the  ^ater  depth  ratio. 
This  atialogy  is  quantitatevly  correct  for  a  fictitious  gas 
hr^viiag  a       -  %/^y  ^^  ^•^^   ^Hydraulic  ju»p^  in  '^sLt^r  flow 
i»  tisod  to  represent  «0o«pre»sion  shock »»  in  air  flow. 

In  ord^r  to  4e»cm^trate  Uw  usefiUjciess  of  the  ,i£  ter 
t&ble  as  applied  to  tfe«  stisdy  of  supers ojiic  air  floif,  it 
was  origiaally  iRtend©^  to  «heek  &  few  of  the  existing 
principles     of  supersonic  air  flom  theory  an  the  is^ater  table. 
Experiments  were  also  proposed  to  iBvefitigmtCf  that  purt  of 
high  velocity  air  flow  .^bm;t  *t?hich  little  is  Imown  at 
present.   Due  to  difficuxti^si  eBtCouiitertd  in  tho  eoiistructiou 
of  the  «ater  table,   little  tiae  wes  available  for  model 
©oris t ruction  snd  test  work*  "Therefore,  only  mm  ar>t^licr.tiofi 
of  zhQ  w^xcr  LEDie  xs  pmi---  This  is  an  inv;  .  ■ioi^ 

of  high  v@loeity  flo^  ^^roimd  a  fifteen  degrtje  ^#age. 

This  work  was  done  durlnr,  the  aorj-ths  of  February 
through  May  1947  t.L  th*?  n«ii^.s»-P^ii'^«?r  Polytechi:d.c  Institute, 
the  mater  tabl®  wss  constructed  i»  the  Mechfiriieal  Eiagineer* 
iiig  Labors tory  of  th<^  Eens&elaer  Polytechnic  Institute. 


-? 


?ABt  I 

THE  J^AtSE  AHALOaY   OF  iA  .  MIC  iOH  FLOS 

BeverslblHt^r,  in  the  eas«  of  &  mem^Tm9siM,e  fluid 

ftii^  as  air,   is  dafizt^d  Wjr  the  st«t«8ieiit  th&t  the  ohazige 

in  iBtera&l  energy  (C^T)   is  G«i«s«d  entirely  W  expansioti 

or  com^resuXon  irork  (IpdV)  or 

J 

#(C^f)-/-P^V  X  0  (1) 

Tte  state  of  »  eoBtpressible  fluid,  is  defin^t  by  the  gas 
equ^ti(»ai 

FT  r  RT  (JS) 

If  the  gee  eqoatic^  (Si)  is  aifferentiated  we  havei 

sji^  since,  at  a^  inetent^  Y  s  M  t   thim 

or 

mT  -  KdT  -  iaS£  (3) 

gabs ti tut  ing  this  value  for  PdV  in  eauatioii  (l)t 


but 


^ 


C  ^  Fj.  C     -     y    K  i^iere      y*  ^j 

S  T-     P  -  "Jul  I                  '     ^ 

equation  (4)  becomes j 

r-1  -^ 


dT^  y-i  a?  (5) 


or 

Z-  -  OL)^  (6) 

the  gas   f?-*ti^t.lrvr:.   u^icg   f^l    -^-otaids  of  alri 


P  r/^gBt  (7) 

0si&i  this  F  in  equftticm  (0)  gives  t 


Hr  eotabiiplttg  e^iU&tiacis.   (0)  &n4  (d)f 

»|li&tlons   (S),    (8),  toid   (9)  define  Uw  rfciatioii^hip* 


b«twe«a  P»/^»  tuui  f  for  coK^ressible  gas  flov  whieh  is 
reversible  (no  he&t  2^dae4  &i^  HiO  frietion). 

'Stm  g«a«ral  energy  eqM&tXon  for  gas  fXowf 

dQ  -f.  dl,  -s  CMT  -f'-lPd?  -hJifa?  -hj^rav 
For  the  Q€ise  of  no  he£:<-t  «^ded  axul  no  outside  work  dotie  on 
tike  g&s   (dlQ  «  0,  dl*  »  0)s 

^^  J  J  Jg 


but 


Hieti; 


CydT  <»  _iPd?  «*.  ivd?  «     C,.dT 


Cpdf  +  ^V€iv  tt  0  (10) 

^  Jg 


Integra tiug  l>*gtweeB  tli«  liiaits  of  T^  sjsd  ft 


If  To  is  taken  at  the  stat©  of  rest     (»<^  r  0),  equfttioa 
(10)  bee<»iesi 

08 lug  Cp  in  the  waits  of  foot- pound »/oefr««  "^JuakiMmt 

IP        2 
Cp(To  -t)  =    i|v  or  ▼  r  ^^gCp(To  -t)  (1^) 

the  energy  ec* nation  In  th©  c&»«  of  frictionles* 
w&ter  flow  states  that  the  s«n  of  th»  potential  energy 
aiid  the  kln&tie  eatiergy  of  ft  water  partlele  1&  constiint, 

€<»i8id6r  i»  flow  filament  (Fig*  1)  wlii<^  ^fii6S«s 
tlirough  ths  point  j^,  z^  of  the  Initial  erosB  seotloe 
X  -  0»  Along  this  fllaiHBit^  l>et«e«&  tli«  pressura  P  «I3^ 
the  ▼elocity  v,   th«  energgr  a  isstion  ist 

P  +  i,or    +^gE  -  Ccoistant  -  P^*  t^v^+pg»^(lS) 
OR  the  fr«e  surf&ee  of  the  water  I*  is  ccmstant  »ndi  «qx2al 
to  the  atftospherio  pressure,  p&.   In  ih*' t  follows,  we  »&y, 
without  fcrror,  set  this  to  sero  since  only  differences 
in  presstire  i^re  of  physieal  signifieano^  in  th#  cas«  of 
incoiEorcssihle  flo^s»  If  th€  rr^ter  floy^s  from  an  infinitely 
*ia«i  i3i;;siii  thtm  Vjj  c  0«  fh6  curv^tur©  of  th€  fret?  surffe© 
Is  also  i5«ro  at  this  point.   It  is  logical  to  ehood«  this 
point  as  the  reference  point  x^,  y^,  z^.  fh«  corresponding 
winter  depth  is  aenoi^eo^  ^  ^o*  ^^  ^^*^  sbove  reference 
pointy  th^  Bamoiilli  ©  i^istion  ist 


or 

^  «  2g(%  -^  if^  -P)/^  (14) 

It  is  reasorLi'..bl«  to  asstawi  that  th^  vertie^.! 
aeeeieratioR  of  the  iii&ter  i*  n^gliglbl^  eompar^d  isitli 

th«  accelrrc^'tion  of  srr:;vity*   nidcr  this   A^-.im-:tian 
tli«  »t*iic:  pr^tsure  at  a  potnx.  m  Liic  lacit:.  ox  i'lom  Varies 
linearly  Viitii  th©  ir«rti©sl  <iistsBe#  from  th&t  point  t6 
the  free  surfaces 

P  r  f g(h~^)  (16) 

Sebstituting  (15)  aiMi   (16)   in  (14)  gXv^&t 

V^  «  3^g(fe^  -fc)   =5  ^iAh  (17) 

tlie  above  energy  eqii^tioB  ts  valid  for  the  tl&m 
ttl^M^mt  passing  throiigh  y     aM  s     at  x  r  ince,  at 

X  r  0>  all  tlie  fllas^Bts  srhich  lie  one  above  tlie  other 
have  the  ease  h     ^ind  v     (zero),  &I5«1   since  eqiiatioii  (17) 
does  not  contain  z,   the  velocity  \  nd  y  is  eof^stamt 

over  the  entire  depth  aiid  is  eqiisl  to  thfe  differ eiace  in 
height  hmtrmmi  the  total  head  (1k>)  ©lid  the  fr^e  level 
(h),    (4h),   ....  i.^.^-.  L  J  -.    «?-,|ttal  h  *  The  aaximaai  velocity^ 


therefore  Is 


V  -  1/ 2L 


If  we  write  eruatioii  (17)   Id  diaensioiiless  foms 

Ah         ,     h 
'»ax  ^o 


<e^>'  -e--  ^-i  ^^^> 


-II. 


Trem  %quBti<m  (IB)  it  is  obvious  th^it  th€  lai&jLlAai 
velocity  i»  j    gas  is « 

r  |/-gC  To 


aad  (^~?~)^  =  t-^  ^1  ^  P  .r    1^  ^  (19) 

FroBi  eqm  tions   (18)   tXia  (IS)  it  any  b©  se«&  thgit 
the  ratio  of  the  velocity  to  the  aaxififfua  velocity  for 

the  WRter  aiid   the  ^ns.   flows  becomes  ^m^l  if 

or 

Etiiiati<H5  (;  0)  shows  that,  tarith  respect  to  ^h^ 
velocity,  there  exists  sii  anslogy  between  the  tv^^o  flop's 
if  th^  depth  ratios  are  ctaapt^red  mith  absolute  gas 
teaper<&tur©  ratios. 

5y  eosparing  th©  © quatioR  of  continuity  for  the 
_Y   watdr  flow  tilth   a  free  surface  to  that  of  &   g&s,  it  is 
possible  to  obtain  another  aiw^logy. 

Consider  at  x^,   a  stt^dl  fluid  priss  »ith  th« 
diasensicms  dx,  dy,  and  h  (Fig*  2).  Let  «,  b^  aiJd  e 
represent  the  ccospcHients  of  the  velocity  in  the  direction 
of  the  X,  y,  and  %   axes  respectively. 

If  th©  ss«imption  is  Bade  thr^t  the  vertici^l  accelera- 
tlcm  of  the  wf^ter  is  negligible  in  comparison  to  the 
acceleration  of  gravity,  eunati<m  (16)  can  be  used  in  dif- 
ferential form  ast 


Tha  right  sid#»  ojT  tfee  i^bov^  ^qwstlmm  ar«  iiKJ#p«xHient  of 
«,  whleh  m^xm  that  the  horl^oi3v...j.  *.ae#ler«tior*s  of  ail 
l^oiats  al<»ig  a  vertical  &r«  mlso  ijid.ep«t^«iit  of  2.  %« 
horijicmtai  e^Map<menta,  »  er^d  b,  are  thaB  cojistant  oir^r  th« 
depth  h. 

fha  continuity  ©qB:atlon  for  this  trp®  of  w&t#r  flow 
states  that  the  rate  of  sulss  flaw  into  ttm  prism  is  «itml 
to  the  TH%e  of  flow  out  of  th^  pris».  Slue®  the  d^uslty  of 
the  »at#r  Is  constant,   the  inflowing  voX^yEse  per  uBlt  ti»« 
(dS^  )  wetMt  «qtial  the  outflo^^lng  voXiaa^  (4Q  ^  Am 

Is 

%,it  -  ^^  -»^*x)  (h*i|dx)^  -«-  (h4^){h4^||dy)4x 

W^  exp&Kalng  ar.d  n#gleetlr,^  Ir finitely  smp..ll  magislttidea 
of  higher  order,  w©  haire$ 

and  divii:?ln^  toy  dxdyj 

Eqiiation  (1:1)   is  th@  co^tlniilty  @4uatl<m  for 
stationary  ^Jrter  flo--* 

the  coj:it.ii:iUii.y  i;.,iiiitloii  lor  &   t*-*o  diffienslon&l 
ec^presslble  gas  flo'#  Isi 


"i  I'iTi'  All      <^       11*11  it  I    wHi/u     o»    fj  s  -'--^■i  i 

Equations  (El)  liBd  '{m)  obviously  hfive  the  9^m» 
forau  Fr^m  th^s©  0  ti^tlons  it  is  |K>s&ibife  to  eeriire  s 
fiy?^©r  condition  *.^.    ^i...  ..x....^v*gf'  of  th-^  ^  .'.    rj^ 

n&meXy,   that  t^e  ol^nsity  ol  th©  gas  flo«  eox-r^sfocids 
to  Ui0  depth,  h,  of  thf;  waiter,  Ex;pri?ssiiig  this  aiialogy 
ms  a  diiseasiojale-ss  ratios 

-t.  r  ^  (m) 

It  is  BOW  posiiifele  to  ii::rtfestig&te  the  ptiysleaX 
aatiir^  of  this  gas   ^h'-^n   .e  ^re  comparixig  with  the  £1  -- 

oT  watar  ^itli  &  frmm  ihurrac#«   The  ^di^^fcatic  aq\aatloii   (§) 
states  thatt 

From  (hq)   cIkI   (£5)   we  hsv^i 

T         h     ,..^^    f  M 

o       ho  C,        lio 

^f  substit'Stiiii;  these  values  In  e-rusition  {B) ,  we  have 

the  e.imcioi.5  ^ 

mhi«h  ohvi(msly  is  satisfied  <miy  if 


iMMNAi|M««M*i^aaHMHv>«MMMMBfwn|i 


Thms  it  is  shofm  that  the  flo*  of  the  water  is  i:uiuBti* 
t&tively  eoMDarc:,ble  with  th<^  flaw  of  a  gas  h&viiag  a  ratio 

mothmr  msiitlog^,  probably  th«  i^ost  iia^port^xit  from 
&a  experia^sttal  point  of  viei^,  1^  ohtiiiji©d  frcai  tha  ga» 
sqij^ticHSf 


or 

1^  ^i^stitutisg  th^  values  of  J^  #r^  _-l  it^  ttiras  of  -A- 

B=qttsttoii  (ii4)   sac  sl:;o  bfe  €i«riire^  by  using  either 
of  th«  adii^betlc  e-B^  tlcn§   fe>   or   {9>   sTsd    JT*  „;.o. 

Par  a  eoi^r^^^^^^uj^'  ^       .    ig^^i^-y  of  propagstioi* 

of  a  soi2»di  wfeTe  r«   -  l/y   '    .  ^  v^elocity  of  |iror>^?^v-' 

tlon  of  SL     surfa.«#  '?^cre  in  ©h-r.Ilow  water  isi  ir|     ~  Y  gh**^, 
Ttiis  velocity  is  e^x/. --..vi  the^cri-.^a^j,  t^^cri-^^city*. 

In  a  gaa  ejt|?iiiiidiiig  tiirot:^h  a  Bosxle  tJcierie  is  a 
critical  pr^ssiiTt  ratio  ??t  »'lil?^h  s<mio  veioelty  o««mrst 
at  the  throi^t*  In  the-  e&se-  oi   -.^^    -,-^w*i.  :.     )^«   '.. ...-w.;-, 
til®  critie&I  p^essui'ts  r?-itio  is    ^  »  l»Bte*;,  If  th*  rdio 
of  the  pr€(&s\xr0  at  fch@  nosxle  efstra^jso^  to  tfe«  thi'ost 
r-re.'«^iire  (fji..^-'   ^  -    '^^-ss  thfe»  tM»  r&Xne,  sonie  velocity 

*    X 

cusm^t  o^aur  ii  ;k>^z1«  throat*  An  atml^gous  c<KB4i* 

ti£»  oeo'iirs  in  tim  waiter  flow, 

ffe©  .->'tt:fr-.  :i   ^i^dociiy  of  «--'3t«r  in  /gh*  Fr<m  th« 


1»     Keil  P.  ^iley^  *Th«  nieraodynj^^slns  of  sir  -t  Hli^ 

¥e loci ties ♦• 
ii*     Ea»uit  F^rei'        "  -tioi.   ji    u^u  ii«t.£««is  of  O^s 

SyasiKUs  tu  ..-^^_.   ^^,--    -^iUi  Fret,-  i.urfae««,  Part  I, 

^^^CA  tH  0S4,   |i*  S,   1S40. 
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Br  substituting  ^  :::   V  f^$   ^qu&ticm  (17)  h^oamm^ 

or 

h     _    E 

Thmx  at  «uEQr  point  where  the  water  d#r>th  is  two  thirds  of  the 
totsl  hesa,   the  water  i«  flowing  critical  velocity. 

In  sir,  if  the  l^ch  nt^jiber  (M  r  ^)  -^^  ^^^^  than  one, 
the  flow  Is  said  to  be  «aubsoiiic»»  aiad  if  W  la  greater  than 
cm©,   the  flow  1»  «mis>ers<^ie»*  t?    ^^t   r,  if  ~  iis  less 
th&u  one,  the  flofe  is  kiioim  as  «&tr&a«ir*g  fla»«  ai»J  if  ^ 
Is  freater  thsn  one,  the  flow  is  i'ahootliig  flow^» 

If  v-.ir  is   floT;lrj?  st   j^urersoni.^  velocity,   fi  55_saonti- 
nuity   ii^,y   occur  in  fehiQii  uiij  Vi;iDCit.y   .         •. tiiy  a^QXtt&i*fes  1.0 
&  subsi^nic  value  Khieh  satisfies  the  conditions  of  flow. 
"Riis  diseoiitinuity  is  c-.lled  s   »:l?.ne  ctH&pression  shock^ 
and   wieor^Lically  tctii«?s  pj.sQ^:  iii  i^  yt&ry  short  diii.tiincfe^.  JuJit 
&»  the  ir^loclty  suddenly  decreases  through  s  plane  shock, 
the  pressure  sm&OMnXf  incr@s»e3«  Whes  supersonic  gas  flow 
is  forced  to  cl^»nfe  its  flow  direotion  due  to  sn  ob.--;,*  v^ucion, 
BXi  »angle  shoek*^  will  ooe«r.^     It  is  knoim  thut  an  analogous 
discontinuity,  called  tlie  «hydr«iulie  Jtiai^w,  »ay  occur  in 
shooting  water  flow.  .Asln  a  6^»,    two  cases  of  the  ^xjaa^  are 
possibles   (*)  In  the  right  ^lagp,  shooting  water  is  con- 


5.     Keii  ?•  Bailey,   ^Abrupt  Energy  Tr&nsforaatior^  in 
Flowing  Gas©s«,   Appendix  B,   p.  2* 

4,    Ibid,  p»  ;^u 


verted  into  str easing  flow  aod   (b)   iia  the  sli^jnt  or  oblicis* 
J«»p,  the  flow  ftey  or  laay  not  go  ta  str^aidug  «f  c«r  th« 
Jii»Py  dep€i5diTig  on  the  — ^  of  the  -*ater, 

A  5UUU  ly  of  the  flo«  analogy  Is  as  followrst 


1.1   HID  FLOr  WITH  FHEE  Sl^^F 

Katia'e  of 
flo^  medium 

IiiCOffipressible  flui^   (wate 

Oe^Bietrlcsjlly 
similar 

..ia^i  bound.iry  vertical 
Bottoffi  horizontal 

AiialogQUS 
sagaitudeg 

.  .  —    ^^_^— J 

Velocity  JL^  ,^ 

max     g 

tempera tute  ratio, «5^ 
Density  ratio^^-I- 
pressure  ratio,  -^ 

Fq 

Doiaj4  velocity,  v| 

Smb sonic  flo« 
Supersonic  flote 

Compression  shoek 

(plane  &r»d  &ngle) 

Velocity  ^  ,^ 

Depth  ratio,  ^ 

Eolith  ratio,  ^^ 

i^fuar^  of  <i€.pth  r&^tio,^-?!--^ 

Wave  velocity,  v| 

?  .Rch  nu«feer,  — i 

ir treating  flow 

thooting  fion 

B^dr^ulie  ixmp 
(norisal  &i»3  slfiJit) 

1-/S 


tkvm  f&r  wm  Ja&v©  siiiown  the  analogy  thnt  can  fe« 

betw4»€^  the  c<»ipr#»^icm  &hoe2c  ©jod  tfe©  liydraulic  Juap, 

feowever,  do^a  not  strietly  ho3-d.   The  eii«?rgy  e^uatlo©   (17) 
|»«tw##a  the  velocity  fcnu  uii^  di^t^iii  i^f  *»4it«r  flo^?  lat 

where  tfe©  tot&I  h#&<^,   (ho   )^  ti>  eonstajat.  lis  thm  cas©  of 

th#  hydr&tilic  Jum^  &  laortior*  oi    ^.ut:  ^ii.^>i.ia  imtsrgy  ar  the 
%Fet«r  Is  »oavert€«i  into  heat*  For  this  reason  the  total 
he&d.  ©ft^r  th€  jiiup,    (h^),   i*  smiller   than  the-  total  head 
h&tQr&  th€  4\a»|>,    (feQ   }•   After  the  jimp  thfe  euer^  eqtmtioa 
1st 

IT     r^2g(h^-h) 
the  mtmtgj  io^s  during  th*s  jump  b#rirs  &  simple  ral^tic^ 
to  the  intensity  of  the  jx^p*  In  the  flow  aver  a  horiz-an-^ 
tal  bottoa  tlie  pot«sntl&l  energy  is  a  aliiifiniM  if  th«  water 
^fepth,  h,   is  i^@ro.   For  a  »ass  of  *rt-er,  m,  at  a  rit^'^th,   h, 
the  pot^^ti^il  energy  i£   F  r  Bgh/^» 

Sisae  the  kinetic  eaergy  at  :poiJat«  of  r^st  is  eqii^l 
to  Eiiro,    the  Allergy  ios#,    (^e),   -^hish  ooeurs  In  the  hydraulic 
jiasp     aay  h«  ec^putea  <-.b   th^  difitirejicc  oi   1.1.4.  pati^iitiax 
energy  at  kk  point  of  ^4iro  velocity  before  ^Osd  after  the 
i\mp,  or 

By  dividirig  this  energy  loss  hy  th«5  energy  before  the 
JiSBtp,  #  s  a^^#    -H«  relative  ©nargy  loss  is  obt^iined  ast 

E  *  Ai  •  1~  hl/ho  (^®) 


This  is  the  relstiire  ajiuuj;. ..  of  efiergy  which  h&s  hi^im 
ooftverteci  into  heat  snd  is  lost  energy  insofar  &s  w&ter 
is  coipicerned* 

Ic  4i&  giLi»y  tho  h«&t  g«£Ler^t#d  by  a  shock  is  i^ot  lost^ 
but  ts  m»T^lLy  ooBverted  into  thtr&al  mn&rgf,  ftad  the  tot&l 
tft»^rature^  coiiiiequejatly  tha  tot&I  ^mergy ,   is  the  &nM^ 
befor«  Djid  «^ter  the  &hock«  Slnee  th#  cuMnparabXe  m&tet 
s&gjaitu(i«  to  the  gss  t€ai»erattir«  is  the  d«pth,   the  analogy 
is  a©t  strictly  true  aftejr  s  hydraulic  ;JiaBp»     The  eiiergy 
loss  in  water^  is  «xtreaely  siai&li  '^v  t-       targe  region  of 
^\wp  intenaity,  the  relative  loss  b^iixig  less  than  onm 
percmtt  for  Jt  up  to  S.O.  As  a  result  of  this  s«all  sho«k 
©aergy  lo&s,   th^      Yi-Iogy  of  the  two  types  of  flov    is  still 
VfiLiid  as  &  very  ciosa  &pproximati<:m  wlthia  the  r^m^ e  of 
Wejsh  maalbeTS  currently  employed* 


5.     Prels«f«rk,  pp»  cit,  F^rt  II,  p*  **il-;^ 
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PART  II 

CO^F^TrnrriCB  OF  TKL  i..  ,I„  ^;AtBK  TABLE 

fhe   essential  element  of  a  water  table  Is  a  horisioxitaX 
surface  vvith  vertical  sides  over  which  «  shalloir  stre&n  of 
w&ter  flows,  the   necessary  adjunct*  to  the  surface  are  a 
freii»«.ork  to  support  it,  the  tanks,  p%mp,   and  piping  to 
suppljR  the  flowing  ^ater,  am  the  necessary  sie&snring  instru- 
«^at».  The  following  description  of  the  water  table  construc- 
ted in  the  Mechanical  Engineering  laboratory  *  ^he 
Rensselaer  Polytechnic  institute  will  serve  to  illustrate 
one  method  of  fabrication* 

Ifie   surface  is  a  sheet  of  r,late  glaas,  i^   x  48"  x   ac*. 
Olass  was  seieci^cG  lor  Li:.o  I'^^iions.  First,  glass  xs  one   of 
the  ffiost  frictionless  surfaces  obtainable  sufficiently  rigid 
to  Balntain  a  horizontal  plane,  fecond,  the  use  of  glass 
aakec  possible  the  taking  of  photograpn^i  oi    \>ii&   flow  either 
fr<»  above  or  unoerneath  the  surface  of  the  table,  illualna- 
tion  being  supplied  on   the  side  opposite  to  the  camera* 

The  frase  is  constructed  of  white  pine.  Tht  tc^ble  box 
is  constructed  of  2"  x  6»  plank,  the  legs  of  4«  x  4«  timbers, 
ar^  the  crossbars  and  leg  braces  of  *c"  x  4«  stock*  An 
Isoetetric  drawing  of  the  fraae  is  shown  in  Fi^*  «  #  ®^ 
actiial  dimensions  can  be  foiifid  in  Fig.  4.  All  joints  are 
aaortised  end  secured  by  bolts*  To  provide  for  leveling  the 
^   table,  the  bottom  of  each  leg  is  fitted  with  a  i»  steel 
^    i»X&te  into  vhich  is  threaded  a  §«  steel  bolt* 


Th&  j$l&»n  is  stipiiort^c  Qu  hhts  sides  hy  1^  x  &;«  strips- 

rest  i»  grooves  cut  in  the  ^^  x  4»  erossbars*  ^t^lls  of  th« 

support  arr^^sgeaaBt  are  ^kuv.i;  xx.  rig.  ?•  !l^o  i»  ^ag^.^;  ironfe, 
surfaced  witte  mmmni.t^,  placed  oae-i;lilrd  Ui^  length  of  th» 
$!$.&&  froii  eaeh  end,  siipply  si3p|K)rt  across  thc^  tabl#,  *rh€: 
glass  is  heXdi  firmly  in  plaec  by  i^  a:  i:'^  v -jv^v   ^i^xxp^  rumxxig 
tfee  leiigtte  of  tht  t&^l©,  th€-se  strips  are  covered  wit^ 
eopjB^r  flasiilTi^  iii  oraer  to  proride  SBOoth  ¥#rti©ai  walls 
&^  to  prevent  the  possibility  of  i^rirping  due  to  Rat^r-to- 
«ocKt  oontact*  To  insure  -^^tertightiie&s  &tid  avoi4  0.  setsl- 
to-glas^  eont«.et,   a  strip  of  Jw  ^i>o»-g#  rti-febisr  is  seoi^red 
fe«t^#€m  the  aopper  &i5d  the  glass* 

*Bi^  t&iaks  are  coiiBtriieted  of  l«i*gfcuge  roii^scl  sheet 
steel.  Dimensions  of  the  forward  tanlt  &re  shovm  in  Fig*  5, 
This  t^tk  h%^  A   f^-t>^elty,  up  to  th^  lip,   ot'  --^fenut  five  cu^ic 
fe<*t.  The  fsrwi-ra  tM^fc  ib  boitea.  to  the  fr^^si^  qb  thre^  ^i<les. 
The  fo\a:'-irich  lip  rests  us>an  th©  forwc^rd  E»  :i  4«  cro&sher 
aiirf   T?rDlc:ats   qvv:'^-'    th^;    fnn.^wd   &ds:e   of  tha   ^l?SiS,    The   IIg  is 

|«  ^pomg^  ruhh4^r»  A  thr#e-i^ch  flange  with  &  staiadard  li« 
t>ii5e  thread  is  brased  to  th©  bottom  of  the  tBnk  to  aacoa^^ 

Diafensians  of  the  after  tank  &m  shoi^ia  in  Fig»  6.* 
It  bB.s  a  e?.p^.citT,  ut-  tc  the  lip,  of  about  5.S  cubie  feet* 
fh@  after  %cv.^-  -.   1.0   wn'.-  ^-^x  i.fc;j-  ^ad  of  the  fraJMi  and 

th0  <m«-ineh  lip  rests  ir*  a  groove  cut  in  th@  ftfter  «i«  x  4" 
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crossbar*  A  stTtp  of  copper  flai^hiBg  rtms  tr&m  ^4®  after 
edg«  of  the  gX&&&p  o^mr  the  taiikli^,  ane  into  tfee  tank*  this 
provides  &  ssiooth  surfiiee  for  i,uss  *;;  i^er  liow  sfict  pxeveutfi 
«rater-to-.®oo4  eoataet.  three  fomr-iisoh  flanges  with  st^Bdard 
two-iBch  pipe  thraad  are  hrased  to  th«  bottom  of  th^  t&nk  to 
aeeomfto^itte  th#  di&eharge  pipeliii«s, 

Tha  water  for  this  table  is  self -coaataisedi  aM  recir« 
ouXating*  For  this  purpos^e  thera  is  &n  80-galios  capacity 
»\aip  teslc  at  the  afttsr  ©na  of  the  table.  Th#  piping  ssd 
connectioiis  b«tw«€sR  the  p\mp^  t%a  table  tasks  ^  and  stmp  tatik 
are  sho%^  in  Fig.  4.   Since  the  pxmp  is  an  alternatii^  ciir- 
rcmt  «l«ctric-driVi;n  coasttiat  spm^  centrlfiigal  tyi«,  a 
ocm trolling  valve  on  the  discharge  side  of  the  ptURp  aM  a 
n&ter  by-p^ss  line  ^re  u£i*d  to  throttle  wad  divert  the 
flow.  All  |ii|5i^  is  two-inch  steel  excr^^^t  for  the  ^-p&sa 
llTie  aBG  short  se^ctioE  leidicg  into  UiC  ror^vara  taaJc,   which 

are  li»  steel*  Valvas  arc*  loc^^ted  in  tlie  three  discharge 
lines  iram  the  Rfter  t^iik  to  regulate  brv.ck  pressure  on  th^ 
tabl^  suri  .  ■.., 

to  measure  Wf^ter  flow,   th^r^j  is  a  flat  plate  orific«^ 
l»fc5»  in  di^isieter.   Id  the  It^  line  les^ding-  to  the  forward 
t4&XUU   i-i'O^surt;   tap!*,   iocai^eu   in  a^sccr^     ^-m::^    -ith  th©  A^iE 
research  publication  on  fliiid  met€rs^,  are  connected  by 
richer  tubing  to  t  l^-ttjb©  r^^^ter -filled,  m^^rometer..  The 

OrifiCfc     t.nC    in;.i.lUint:  L^.-r  C^-.XJ.l>X'>:.  vtry      Uj      i.v-xghing     th€ 


6,  »FluiQ  M&ter&t   Thair  Theory  and  Applicati<m,«  Part  I, 
ASME  Riitsc^^.TCh  publication,  19S1# 


*»oyiit  of  water  flowing  jua  a  measured  InterTftl  of  time. 
A  caiibratlcm  ciJrve  is  sfeovm  in  Fig*   S» 

To  ®«jr.sure  depth  on  the  water  teble  s\irfaee,   &  hook 
gau^^,  capable  of  ae^surisg  to  X/lOQQ  foot  is  u»ea»  The 
hook  was  r^pXgicm  by  f:   str.^^ight  pointed  rod  for  greater 
conv^nitsnce  anxl  aceuracy*  The  gauge  is  isount'id  on  a  woodea 
block  iBhich  slices  fre«sly  aloijg  a  g«  x  g»  wooden  stringer* 
The  stringer  rests  on  the  table  frame^   Depth  In  me&9wr%d  Iqf 
t&kiBg   the  difference  of  &  readliig  at  the  glass  and  on«  at 
the  water  surface. 

A  wire  b;asket  filled  f^ith  sBtall  roeks  is  lu$ed  to 
ssiooth  out  the  water  flow*   This  basket  is  fsbric&ted  of 
two  strips  of  i^  wire  mesh  ^hich  are  three  inehes  apart. 
The  basket  rests  in   tbi>  ct-nter  of  the  forward  tank  and 
ruiis  trui  HiQZh  or  tnc  t^im.   It  also  serves  the  purpose  of 
&  water  filtor*  Two  iseirs  are  provided  to  su^j>ly  th«  neces- 
sary total  hs?:d.    They   are  fabrics* te-d    of  14  EBum  sheet 
brass,   roliea  fXi-.,    bent   to  xhe  aesirt-     -  r.-  ■, :  .    .i...     unacriiaath 
side  of  each  weir  is  fittod  with  three  e^jually  spaeed 
wooden  blocks  into  ©sch  of  ?hich  j  a  J»  bolt,  the 

woodfe-n  tJiOL'iii   ar^   at^j-u  in  pxace  by  s»a.ll  brass  je& chine  bolts 
•tshich  are  couritfcrstmk  into  the  top  of  the  i^eir.   the  bolt 
heiids  itre  cover  eta  with  soldc^r  ?hich  rtss  s«ooth@d  to  eonforv 
to  the  shejie  of  the  suriac:   of  tht   v;«fir,  The  seouriBg  bolts 
pro^^ct  through  holes   in  the  lip  of  th©  forwv^ir^  ta&k  «n& 
crossbar  and  are  secur^}d  by  nuts  b^i^^ring  agai&et  the  under- 
side of  the  erossl»ar«   water  tightness  is  naintaieed  b^r  iielBi 
eeixlking  wm^mmd  betveeB  the  weir  and  the  t£.ble  sidee  afid 


The  woocifeG  models  used  thus  far  1b  canjimctioa  «lt& 
•a(.|?€riji«at*  cm  the  t«.bie  were  fsbrle&ted  of  i^hlte  pija*  &&£ 
finlslied  with  si?iir  varnish^  ^^tud^d  to  a  fine  polish* 

To  miiiLmlzi    surlacc  tension,   &  r^etting  eg^^t  is  usiid 
in  this  «raLter. 

Fig*  10  is  ^.  photogrypfe  of  the  completed   table  in 
vhich  th«  i'v$.WB,  pximp,   piping,  v^lv^iG,  @£a::oKieter,  an<5  ^epth 
gauge  are  Yislble» 

Fi^.   II  is   :-:   t>hotogra:?h  nf  the  t^-^hlr   top  lookiRf  for- 
wara,   j_xi  v^iiicii  vn-  ci,    ^Ltr,   ana  iroeiw  inter  art*  visible. 

A  typical  experiasBtal  sat^up  is  shown  in  pi&ee  on  th* 
gl&ss« 

Xxi  uiii  coxLTse  oi   construction  t^o  B&jor  Giriicudes 
««r«  encouat«jr«d,  mounting  tha  glass  K»d  obt&ining  laJBi»&r 
flow, 

Tht:  ori.j^i.r>-.x  -uii.^xi^j.oii  v. !i£.   to  u^c  i^  plate  glass  hath 
for  rigidity  nhd   strt^ngth*  This   glass  could  not  be  obtained 
feithin  the  tijae  available.    Pft^.r  inTf^^tir^citioa  of  strt^ngth 
and  rigidity  properties,  ^  dateriRiaM  fea 

be  satisfactory  and  we.s  instttlieo.   This  gl&ss  was  Buoport^ 
caaly  on  the  tour  sides.   About  one  month  after  installati<m^ 
this  gls.ss  fraotia<^u*   Apparently,    the  glasfe  h&d  saggcii 
slightly  in  the  center  and  the  stress  set  up,  in  tise,  re- 
sulted in  fraeture.  A  second  i"  plate  glass  was  instsiiitd 
and  the  two  aagle  iron  braces  %ere  ad4ed  as  additional 
support.   After  about  three  %e@ks,  a  »»&il  fr^iCture  »as 
experienced  at  tht'  after  corner  of  the  table,  tho  cause  of 


this  seconc   hr^tM  is  m>z  del inltely  Ipsowb.  Mo  imp^Qt  spot 
is  discernible  and  the  glass  is  5>#rfactly  ieval.  This  break 
does  cot  affect  the  oner?.vtiDn  of  the  trsble.   k  sh^et  of  t» 
plate  glass  ka^s   uevn  oraert^a  ixiia  wxxi  be  liistelXea  as   soon 
as  del iverM* 

The  Eerfy^olds  RiJimfeer  of  the  flow  an  the  tpfcie  is  renf 
near  the  critic#-l*   In  or<!er  to  ofetaija  ii.su.rvs.x-  liow,  the 
water  leust  proeeeti  froia  its  entraiaee  i^to  the  forward  tai^^ 
throy.gh  the  rock  filt^^r,  over  the  w^»ir  aM  oa  to  the  tabi« 
with  gradiial  tm'ns.  The  i^eir  surfae€  i^st  be  extremely 
SBOOth  since  ai^  protuher&RUe  will  s#t  up  &  disturfesBce 
and  force  t\m  mhter  into  turfeisl^^t  flow,  the  first  weir  us#d 
«&s  coaistrEtctod  of  wood,  A  really  s»ooth  suriao©  coula  Bot 
he  ohtal»#d  on  tlm  »ood  tv&ilahle,   am  this  iras  di&eiird#4 
In  fairor  of  th^  weir  const  rue  t#c;  of  hr&s»«  L^^inar  flos^;  i« 
ROW  ohtaia«!<i  thi'oiighout  "•'•-    r^r^ctlefel  r^mge  of  flo'^-  'c--  i^-***!-- 
ties. 


FMMS   XXI 

APPhlCATlQMB   OF  TEE  WATEii  t^LE 

The  applications  of  the  vater  table  ija  eozmectloii 
with  the  study  of  air  flow  am   s&ny  in  niistber.  the   water 
table  Biay  be  used  to  sttidy  any  type  of  two-dimensional 
air  flow  which  does  not  involire  a  change  in  the  total 
«®ergy.  Some  of  the  more  important  applications,  in  the 
opinoin  of  the  writers,  are  as  follows s 

(1)  the  inTestlgation  of  the  flow  of  supersonic  air 
in  and  eround  airfoils,  inlets,  nozzles,   and  fliff users, 
Models  aust  be  designed  ht^   fabric* tfcd  for  air  flow  with 

a     y«  g,0, 

(::')    confirciPtion  of  the  validity  of  assusptlons  made 
in  t^hfe  aevelopment  oi  air  xi'.;     ory.  Using  )(m     -...  in 
the  e^iUations,  the  theory  can  be  checked  with  experimental 
results  obtained  from  thp  :^cter  table.  The  wstcr  t^ble 
offers  an  excelicnv  ^  ^  i  ior  SLuaying  the  ixi^Lxe  snomi 
trans sonic  region. 

(3)  as  a  laboratory  apparatus  for  deaionstrating  the 
accepted  Liieory  qi   Vc-rious  types  of  gas  flow, 

(me   example  of  the  applies tions  of  the  w^ter  table 
is  demonstrated  in  the  case  of  flow  of  supersonic  sir 
around  a  simple  wedge«  The  purpose  of  the  experiment  was 
to  check  the  water  analogy  with  the  existing  theory  on 
angle  shocks."^ 


?•  Bailey,  0£.  cit,  p.  4. 


A  9k»tch  of  tht^  «i«dg«  «ui«d  In  the  experls^at  Ifi  shown  in 
fig.   b.   Experimental  data  cslb  be  fotind  in  Table  !•  Table  II 
contains  the  results  calci^i^Lea  ironi  tht.   .  .«    criiseniiax  aata. 
TAble  III  contains  the  mstilts  calculated  for  air  with  a 
y*  2.0.  Sample  calctilations  can  b©  fovuad  in  appendix  A. 
A  comparison  of  Tables  II  rnd  III,   shows  the  eorrelation 
b«twe^i  the  ^ater  table  analogy  and  air  flow  is  within  the 
limits  of  experiaental  accuracy. 


tli«  thftory  of  the  analc^gr  of  t-"o  dleettsicmal  air  flow 
to  the  riow  of  -^/xtc^r  feiUi  a  fr««i  sui^x.  '■, 

rh«  details  of  ccmstruetion  of  the  v&ter  tftl^le,  i 

utilises  the  «&t#r  aaald^  li»  tlm  »tvs^  of  air  flow,  1  > 

13ie  operatlffl^  ^xad  9m  ftp^^Uaatioa  of  the  £&ter  tabl« 
h».vm  hm^^  d«KQiistrAt«d« 

fli«  water  table  Qould  be  esplo^retl  to  a  c^cld^  ad«- 
▼axLtaf«  hy  the  orgAUs&tiond  ¥»hich  tuse  th«  fe^i  $up€r motile 
«ii2d  tiam^ls  in  exl4ita&e«  in  this  country  to  s%vt&j  the 
flow  of  high  Teloeity  sir*  Mueh  of  the  feori£  done  at  pr#si»t 
in  the  vixsd  tuimels  could  he  aceorn;  llshe4  ^i  &  w^ter  table 
with  the  result&nt  eofisl^er«^hle  t^tntXug  1b  tiate  &m  exp^tee, 

fhm  9m^^Bls  in  this  i^ork  h&s  he<^  o&  sui>ersonic  elr 
flour ^  how^rTerj^  the  iirater  table  c&a  be  uimd  for  s^hi^onic 
sir  flote  investigations •  It  is  iH  fine  liboratoir  afpars  ttis 
for  demoiistrEtliig  hydraBlie  sisd  hXr  fla^  theory  trlstiplly. 

Tfa©  oper&tio&  of  the  pr^aesiii  .  c/ucr   usbie  f  oi^c  be 
^[ihancect  ^ith  th€  uee  of  a  dept^  gauge  oapable  of  i6«^asurixkg 
to  1/1000  inch.   A  further  refiu^wmt  is  control  of  ^^^  ter 
flo*  is  n&e^sst.ry  „   It   is  sta^^eiiced  tll^t  the  preseot  ymuy 
end  Kotor  be  replaced  by  e  «©»troll;sble  direct  current 
electric  drive  variable  sfH^ed  eentrifugal  pt»p.   It  is   al&o 
sugsesteo   th.!^t  lxior<i*^r   to  secure tely  sea  sure  shock  fugles 
on  photographs y  &     parallel  light  eyetes  be  oonetruoted* 
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!•     Belley,  Hell  p.  -  »1^  ^^»rm>^^jmm9^cs^  of  Air  at  Hlgli 

Velocities^,  liolvsm  II,  m.  ^,  ^m^m^-  9f.  ^  mmr 

2,     Bailey,  Keil  P.   -  "Abrupt  "Bmrnrgy  TrEHsfonaatioii  Is 

Flowing  G&ses*»» 
3«     BaJklUBet^nr,  B.  A«  -  oHftSracilics  of  Opmt  Cb»mml%m^ 

Wbs9  WBsQr^w-'Blll  Book  CkKBip«tiiy,  K<i{^«  York,   1&54« 

4,  «Fluid  ir^ter,  Tb«ir  Theory  asd  APFlieiitio»»,  part  !♦ 
^err^rt  of  Assoei'tion  of  rech  rIc  1  !>igine©rs  Special 
Qoaa&ittaa  cm  FliiicI  liet^rs«  Aij^  ^es^- ^rcli  publicfe.tlaii» 
y«^  YorJc,   Ij^SI, 

5.  Preiswerk,  g&rnst  ~  wApplie^^tian  of  the  Vethod^  of 
0^»  I^rn^aics  to  ■'  Flo^^  *ai-  ui-i^ioe*'* 

Idrisory  C09SBttitt©#  for   ieror.fiuties. 
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A  m  cross  £«ctloe»I  ar#«  of  «.c^ter  flow,   equal  to  fc«t>   ftf 
M^  s  wtA  Raal^r  of  wst«r  flo»>  ▼^vj 

vj  r  ireiocity  of  propog&Xiam  of  &  surf®e«  wmw,  fpm 
w  m  «i<ith  of  »at©r  flow,  ft« 
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Abrupt  Energy  Transformation 
in  Flo\ving  Gases 


Paper  No. 
46— A-71 


By  N.  p.  bailey,!  tROY,  N.  Y. 


After  setting  up  the  basic  equations  of  energy,  flow,  and 
acceleration,  this  paper  compares  theoretical  and  actual 
plane  compression  shocks  in  tubes.  By  making  the  single- 
plane-angle  shock  a  special  case  of  a  plane  shock,  curves 
for  solving  numerical  cases  are  presented.  Even  though 
most  angle  shocks  encountered  in  engineering  are  three 
dimensional,  it  is  felt  that  the  simple  theory  of  single- 
plane  shocks  is  helpful  in  understanding  such  phenomena 
without  the  complications  of  three-dimensional  theory. 
Traverses  of  three-dimensional  angle  shocks  in  nozzle  and 
orifice  discharges  are  presented  and  discussed,  and  tem- 
perature traverses  through  shock  diamonds  in  high- 
temperature  streams  are  shown.  Although  no  complete 
explanation  exists  for  thermocouple  readings  in  excess  of 
the  total  gas  temperature  in  such  high -temperature 
shocks,  their  existence  is  established.  As  a  last  case  of 
abrupt  energy  transformation  in  gas  streams,  tests  of 
combustion  with  flow  in  a  constant-area  tube  are  given. 
Experiment  and  analysis  agree  that  the  assumptions  of 
constant-area  flow  and  steady  or  continuous  flow  do  not 
appear  to  be  simultaneously  tenable  when  combustion  is 
present. 

Nomenclature 
The  following  nomenclature  is  used  in  the  paper: 

A  =  area,  sq  ft 

Cp  =  specific  heat  at  constant  pressure,  Btu  per  lb  per  deg  F 
C„  =  specific  heat  at  constant  volume,  Btu  per  lb  per  deg  F 
Ci  =   constant  of  integration 
Ci  =   constants  of  integration 

d  =  total  differential 

d  =  partial  differential 

F  =  friction  force,  lb  per  ft  of  length 

g  =  acceleration  of  gravity;  32.2  fps  per  sec 

H  =  rate  of  heat  release,  Btu  per  sec  per  cu  ft 

/  =  778.26  ft-lb  per  Btu 

L  =  mechanical  work,  Btu  per  lb 
M  =  Mach  number 

P  =  pressure,  psfa 

Q  =  heat  added,  Btu  per  lb 

R  =  gas  constant;  53.3  for  air 

T  =  temperature,  deg  Rankine  (R) 

t  =  time,  sec 

V  =  specific  volume  cu  ft  per  lb 

V  =  velocity,  fps 

W  =  weight  flow,  lb  per  sec 
X  =  distance,  ft 
a  =  deflection  angle,  deg 
7  =  ratio    of    gas    specific    heats; 


1.395    for    cold    air 


1  Head,  Mechanical  Engineering  Department,  Rensselaer  Poly- 
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6  =  shock  angle,  deg 

p  =  mass  density,  slugs  per  cu  ft 

Introduction 

Mechanical-engineering  literature  is  rich  in  analytical  and 
experimental  information  on  the  continuous,  steady  flow  of  gases 
but  there  is  need  for  a  wider  familiarity  with  gas  flow  when 
such  discontinuities  as  compression  shocks  and  combustion  with 
flow  are  present.  The  aim  of  this  paper  is  to  present  some 
analytical  and  experimental  information  on  such  discontinuous 
and  intermittent  phenomena. 

Flow  Momentum  and  Energy  Concepts 

In  any  case  of  flow,  such  as  illustrated  in  Fig.  1,  the  difference 
in  the  rate  at  which  mass  enters  and  leaves  the  volume,  Adx 
must  be  balanced  instantaneously  by  the  rate  of  storage,   or 


b{pvA)  4j   ^P 

ax  =  — Adx  — 

bx  bt 


[1] 


At  any  instant  the  mass  of  fluid  in  the  volume  {Adx)  is  (pAdx) 
and,  if  it  receives  a  change  in  velocity 


dw  bv 

dv  =  —  dx  -\ dt. 

bx  bt 


[2] 


the  acceleration  is 


[3] 


dv        bv  dx        bv  dt       vbv       bv 
dt        bx  dt         dtdt        bx        bt 

For  a  wall  friction  or  obstruction  force,  Fdx,  a  summation  of 
forces  in  the  flow  direction  gives 


AbP  , 

dx  =  Fdx  —  pAdx 

bx 


vbv       bv 
bx        bt 


=  0. 


[4] 


A  useful  relationship  may  be  had  by  multiplying  Equation 
[1]  by  V  and  adding  it  to  Equation  (4)  to  give 


AbP   ,   vb[pAv]  ,     ,v  t)t; 

— f-  — =  —{pvA) pA 

bx  bx  bx 


fe)- 


bt 


rpbA~\ 

By  adding     to  both  sides  of  Equation  [5]  and  remember- 

L  bx  J 

bA 
ing  that  —  =  0  for  a  rigid  passage.  Equation  [5]  may  be  written 
bt 

as 

o^  .        r.^  ^[P^   +  pAv']  b(pAv) 

PdA  —Fdx  = dx  -\ -—  dx [6] 


dx 


bt 


b{pAv) 
For  the  usual  steady-flow  case  where  -— —  =  0,  Equation 

[6]  says  that  the  net  wall  force  on  the  gas  in  the  direction  of  flow 
is 

[71 


PdA  —  Fdx  =  d  [PA  +  pAv^]. 


1-^1 


Discussion  on  this  paper  will  be  accepted  until  February  3,  1947 
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where  [PA  +  pAv^  is  the  total  momentum  per  second  in  pounds 
passing  any  section. 

Since     „  =  — —-     and  the  Mach  number     M  =  — , 

r       3«^J  L  VygRTJ 

it  is  often  convenient  to  express  Equation  [7] . 


Net  wall  reaction  =  PdA  —  Fdx  =  d  [PA{1  +  yM^)] .  .  [8] 

This  statement  is  true  only  so  long  as  the  flow  is  steady  and 
it  does  not  hold  when  the  mass  flow  (pAv)  changes  with  time. 
However,  since  any  such  flow  variation  must  go  through  a  re- 
peated cycle,  the  average  value  of  J"  d{pAv)  must  be  zero  and 
Equation  [7]  may  be  used  to  evaluate  the  average  wall  force  even 
when  the  flow  is  unsteady  or  intermittent.  This  point  will  be 
covered  by  experimental  work  later. 

As  gas  flows  along  a  channel  as  in  Fig.  1,  any  thermal  energy 


R 
P 

ir 
V 


-dx 


Fig.  1 


R  +•  d  R 
Pt  d  P 

VtdV 


{dQ)  that  is  released  or  added,  and  any  work  done  on  it  {dL) ,  as 
in  a  compressor  impeller,  must  be  accounted  for.  Some  of  it 
will  go  to  change  the  internal  energy  (CvdT)  and  part  to  do  flow 


if) 


work  (  ^^^  )  as  the  gas  flows  through  a  pressure  change.     With 

a  pressure  change  (dP) ,  a  compressible  gas  will  undergo  a  volume 

/PdV\ 
change   (dV),   and  do  an   amount  of  expansion  work  (   — —  1 

and  an  increase  in  velocity  would  induce  a  kinetic-energy  change 


U**) 


This  results  in  an  energy-balance  equation 

dL  +  dQ  =  C,dT  +  \  VdP  +  I  PdV  +  ^  vdv 
J  J  Jg 

A  reversible  process  is  characterized  by  dQ  =  0  and 

1 


C^dT  +  -  PdV  =  0 


[10] 


Such  a  process  is  thus  defined  as  one  where  the  only  internal 
energy  change  is  the  inevitable  one  resulting  from  compression 
or  expansion  work. 

A  aseful  form  of  the  energy  equation  for  general  applications  is 


dL   +  dQ   = 


{y 


—   -   dl    +   —  vdv. 
1)     J  Jg 


[11] 


For  the  ca.se  of  constant  total  energy  {dL  =  0  and  dQ  =  0),  if 
Equation  [11]  is  integrated  between  any  velocity  v,  and  corre- 
sponding static  temperature  T,  and  a  final  velocity  of  zero  where 
the  static  temperature  is  the  same  as  the  total  temperature  To, 
the  result  is 


ygRT  +  ^-^^—-^  «»  =  ygRTo 


[12] 


When  this  is  combined  with  the  Mach  number  definition 

M-  =  [13] 

ygRT  ^  ^^ 

The    result    is 

^                             ^'o  (total)  .,_,, 

2^  (static)  =  Z [14] 

2 

For  the  case  of  steady  flow 

W  =  pgAv [15] 

and  this  combined  with  the  gas  equation  and  with  Equations 
[13]  and  [14]  gives 

wVTo       „     1^         (7  —  1)        1 

Equation  [16]  is  very  useful  in  all  cases  of  steady  flow  at  con- 
stant total  energy. 

Plane  Shocks  ' 

For  the  simple  case  of  steady  flow  at  constant  area  and  con- 
stant total  energy  with  no  appreciable  wall  friction,  Equation  [8] 
becomes 

P(l  +  7M2)  =  Pi(l  -I-  7M12)  =  Pad  +  tMj^)  =  const.  .  [17] 

and  Equation  [16]  may  be  written 

I  (7  —  1)  /         (7  —  1) 

PM  Jl  +  '-^ — -  M'  =  PMi  Jl  +  -^ — ~  Mi^ 

=  P2M2  a/i  +       "~      M.,2  =  const [18] 

When  (P2/P1)  is  eliminated  between  Equations  [17]  and  [18], 
one  obvious  solution  for  the  equation  of  constant  area  and  con- 
stant total  energy  flow  without  friction  is  for  (Mi  =  M2),  but 
when  this  root  (M2  —  Mi)  is  factored  out,  the  result  is 

i  +  (riil)Mi^ 

M.^ [19] 

■Kir  2     (y-~^) 

and 

27M12— (7  — 1) 

P2/Pi=    "^     .       ',, ' [20] 

(7  +  1) 

This  means  that  for  any  superacoustic  Mach  number  Mi, 
there  is  a  subacoustic  Mach  number  M2,  which  satifies  the  condi- 
tions of  flow.  This  defines  a  plane  compression  shock  which 
is  a  discontinuity,  occurring  theoretically  in  an  extremely  short 
distance.  Whether  or  not  it  will  occur  in  any  case  depends 
upon  operation  conditions.  Since  channel  friction*  causes 
superacoustic  flow  to  approach  the  acoustic,  for  any  initial 
Mach  number  there  is  a  maximum  flow  distance  for  which  shock- 
less  flow  is  possible. 

A  shock  may  also  be  induced  in  a  shorter  channel  by  applying 
the  correct  back  pressure.  For  any  initial  Mach  number  and 
channel  friction,  each  position  of  the  shock  uniquely  determines 
one  back  pressure.  This  is  illustrated  in  Fig.  2,  where  A  —  1  —  B 
represents  flow  that  is  all  superacoustic.     If  the  back  pressurfe 

2  "The  Thermodynamics  of  Air  at  High  Velocities,"  by  N.  P. 
Bailey,  Journal  of  the  Aeronautical  Sciences,  vol.  11,  July,  1944,  p. 
227. 
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y^                                       "^ 

X. 

Fig.  2 


is  raised  somewhat,  the  flow  adjusts  to  it  by  a  plane  compres- 
sion shock  at  some  point  1,  going  to  subacoustic  flow  at  2,  and 
discharging  subacoustically  at  C.     Equation  [19]  indicates  that 


the  Mach  number  after  shock  approaches 


V 


(t-1) 
2t 


as  Ml  becomes 


very  large,  and  it  cannot  go  below  that  value. 

The  magnitudes  of  the  pressure  rise  and  Mach  number  change 
of  actual  compression  shocks  check  welP  with  theorertical  values 
from  Equations  [19]  and  [20] .  However,  instead  of  being  com- 
pleted in  a  negligible  distance,  the  pressure  rise  extends  over  a 
distance  of  from  one  to  five  pipe  diameters. 

This  is  illustrated  by  the  test  data  shown  in  Fig.  3.  Both  wall 
and  center  pressure  traverses  were  taken  for  the  conditions  shown 


and  after  a  series  of  angle  shocks  in  the  first  1.15  in.  of  the  tube, 
the  wall  and  center  pressures  became  the  same,  and  a  plane  com- 
pression shock  was  initiated.  The  Mach  number  as  calculated 
from  Equation  [16]  for  an  initial  pressure  of  14.33  in.  of  mercury 
was  Ml  =  1.37.  Similarly  at  the  final  pressure  of  29.50  in., 
M2  was  0.75.  From  Equation  [19],  a  shock,  started  at  Mi  = 
1.37,  should  have  ended  at  M2  =  0.748,  and  from  Equation  [20] 
the  pressure  should  be  P2  =  29.50. 

The  only  discrepancy  lies  in  the  fact  that  the  pressure  rise  is 
distributed  over  a  length  of  1.3  in.  instead  of  being  an  abrupt 
change  as  illustrated  in  Fig.  2.  The  explanation  of  this  lies  in 
the  effect  of  boundary  layer  and  the  use  of  a  manometer  to 
measure  pressure.  Since  the  flow  in  the  boundary  layer  is  sub- 
acoustic,  no  compression  shock  and  corresponding  pressure  rise 
can  occur  in  it.  This  means  that  at  any  instant  the  large  pres- 
sure rise  across  the  shock  front  is  short-circuited  through  the 
boundary  layer.  This  immediately  causes  the  shock  to  col- 
lapse and  form  at  some  other  point,  only  to  collapse  again. 
Since  the  manometer  reading  is  a  time  average  of  the  pressure 
at  any  point,  a  value  of  one  half  the  pressure  rise,  occurring  at 
1.4  in.,  Fig.  3,  merely  means  that  the  compression  shock  is 
upstream  from  that  point  one  half  the  time  and  downstream 
from  there  the  other  half  of  the  time.  The  shock  is  constantly 
dancing  back  and  forth  in  the  tube,  never  being  further  upstream 
than  the  1.17-in.  point  and  never  further  downstream  than  2.5 
in. 

This  explanation  is  borne  out  by  high-speed  photography,  as 
well  as  by  the  fact  that  an  open-end  impact  tube  tunes  to  a  shrill 
whistle  in  such  a  shock  region.  Consequently,  the  simple  pic- 
ture of  Fig.  2  will  have  to  be  modified  to  show  the  compression 
shock  (1-2)  as  sweeping  back  and  forth  in  the  tube  through  a 
distance  that  varies  with  boundary-layer  thickness,  tube  size, 
and  initial  Mach  number. 

Angle  Shocks 

When  supersonic  gas  flow  encounters  an  obstruction  which 
calls  for  it  to  change  its  flow  direction,  an  angle  shock  will  occur. 
This  is  because  the  approach  velocity  is  so  high  that  the  first 
molecules  to  find  themselves  in  trouble  are  not  capable  of  sending 
a  distress  signal  upstream.  This  means  that  each  successive  layer 
of  gas  molecules  must  run  into  the  same  trouble  with  no  warning. 

This  is  illustrated  in  Fig.  4  for  the  simple  two-dimensional 
case  where  supersonic  parallel  flow  at  a  velocity  Vi,  must  turn 
through  an  angle  a,  and  flow  parallel  in  a  new  direction.  This 
entire  change  takes  place  in  a  shock  front  at  an  angle  d,  from  the 
original  flow  direction.  Such  an  angle  shock  may  be  made  into  a 
special  case  of  the  plane  shock,  discussed  previously,  by  viewing 
it  as  a  normal  velocity  (i^i  sin  6),  going  through  a  plane  shock  to  a 
velocity  Vt  ,  and  this  all  occurring  in  a  flow  field  having  a  uniform 
velocity  {vx  cos  e)  parallel  to  the  shock  front.     This  can  be[|done 


FlNRL  Flow  DiRECnOM 


Fig.  3     Compression  Shock  in  a  Tube 


Fig.  4 
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because  there  is  a  pressure  rise  normal  to  the  shock  front  but 
no  pressure  gradient  to  give  a  velocity  change  parallel  to  the 
shock  front. 

From  Fig.  4  the  condition  that  the  normal  component  of  the 
velocity  must  be  reduced  to  a  value  t;>',  that  will  cause  the  final 
velocity  Vo,  to  be  turned  through  an  angle  a,  may  be  stated  as 
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and  from  Equation  [12]  the  total  temperature  Toi'  is  i 

[27t 


r„/  =  T02'  =  r/  +  ^-^r-^  V,' 


2ygR 
This  gives  the  total  temperature  To-i   as 


=  tan  {e  —  a) 


V\  cos  0 
From  Equations  [13]  and  [14] 

■\/ygRToi  M,  COS0 
1  +  i2^)  M,^' 


Vi  COS  6  = 


V 


[21] 


[22] 


T    '    — 


1  +  (X^J)  M.3 


(t  —  1)  T'oi  M .^sin^e 
+  '-^ ,  =i  . .  [28] 


sIi  +  ^J^I^iIIm.^ 


Using  Equations  [24]  and  [28!  in 

VygRToi'  M2 


Also,  since  t)/  =  i^i  sin  6,  for  the  same  static  temperature  7'i  =  Ti'> 
Equation  [13]  says  that 


V 


Mi'  =  Ml  sin  e. 
From  Equation  [19]  then 


[23] 


gives 


V2    = 


VygRToi 


1  +  ^XZlI)  M,'2 


(7  —  1) 
1  +  — Mi^  sm2  e 


I 

[29] 


1  +  ^ ^  Mi'2        1  +  — Mi^sin^e 


V 


1  +  ^yjiz^  Mi^ 


M/2    = 


7Mi'= 


7—1 


yl/liHm'^e- 


(7  — 1) 


•  .  [24] 


Before  {V2')  can  be  evaluated  from  Equations  [22]  and  [24],  the 
total  temperature  I'oj'  =  Toi'  must  be  expressed.  The  static 
temperature  Ti'  =  Ti  may  be  evaluated  from  Equation  [14]  as 


(7  +  1) 
2  •  2 

From  Equations  [30]  and  [22]  in  [21] 

(7  — 1) 

1  +— 

2 
tan  (e  —  a)  =4 


[30] 


Ml  sin  0 


Ml 2  sin^  d 


(7  +  1)  M^sin  2  d 


[31] 


Ti  =  T,'  = 


Since,  from  Equation  [22] 


Vi   =  Vi  sm  6  = 


Toi 


1  +  ^-^^Mi^ 


\/ygRToi  Ml  sin  6 


l  +  ^^^M, 


V 


[25] 


[26] 


Equation  [31]  defines  the  shock  angle  6  for  gas  flow  at  an  initial 
Mach  number  Mi,  being  deflected  through  an  angle  a. 

Using  Ml'  =  Ml  sin  d  in  Equation  [20]  gives  the  static-pres- 
sure-ratio rise  through  the  single-plane-angle  shock  as 


P2  _  2  7Mi2sin2e  _  (7— 1) 
Pi  ~       (7  +  I)  (7  +  I)' 

The  resultant  final  velocity  v^  is  given  by 

j,2  =  ^yv2'^'  +  v^^cos^e.... 


[32] 


[33] 


P'iG.  5     SnorK  Angle  0  Versus  Deflection  Angle  a 
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Ubing  Equation   [21]  for  v/  and  Equation   [22]  for  Vi  cos  0  in 
Equation  [33] 


^2  =  — 

V 


\/ygRToi  Ml  cos  6 


1  +  — -Mi^cos{e  —  a) 


[34] 


From  Equations  [12]  and  [13] 


M,= 


^ 


ygRTo2        (7  —  1) 


[35] 


Since  vi  and  Vi  are  total  velocities  and  a  compression  shock  is  at 
constant  total  energy,  2'iji  =  2^02 


.[36] 


M2  =  = 


1  +  (x:^  M.. 


cos^  {d  —  a) 


(7-I) 


Mi2  cos^  d 


This  gives  the  final  Mach  number  M2,  after  an  angle  shock  from 
an  initial  value  of  Mi,  when  single-plane  flow  is  deflected  through 
I  an  angle  a. 

Fig.  5  shows  values  of  shock  angle  6,  plotted  against  the  deflec- 
!  tion  angle  a,  for  various  values  of  initial  Mach  number  as  calcu- 
lated from  Equation  [31].  For  each  value  of  Mi  there  is  a 
i  maximum  deflection  angle  for  which  angle  shock  conditions  can 
i  be  satisfied  and  for  any  increased  deflection  angle  a  plane  shock 
I  will  result.  This  maximum  point  represents  the  conditions 
I  where  a  wedge  in  a  free  stream  would  cease  to  have  angle  shocks 
from  its  nose  and  would  have  a  bow  shock  out  in  front  of  it. 

Fig.  6  shows  corresponding  values  of  shock-pressure  ratio  for 
var,''  .y  conditions,  and  Fig.  7  shows  the  final  Mach  number  after 
an  angle  shock.  Whereas  the  flow  after  a  plane  shock  is  always 
subacoustic,  it  may  be  either  above  or  below  sonic  in  the  case  of 
angle  shocks. 

When  supersonic  air  flows  past  a  wedge,  the  angle  shocks  are 
stably  anchored  at  the  nose  of  the  wedge,  but  when  a  boundary 


Fig.  7     Mach  Number  After  Shock  Versus  Deflection  Angle  a 


Fig.  6     Shock-Pressure  Ratio  Versus  Deflection  Angle  a 

layer  is  present,  th  >  shock  pressure  rise  is  distributed  over  an 
appreciable  distance  as  in  the  case  "of  a  plane  shock. 

In  Fig.  3,  where  the  first  angle  shock  originated  at  the  wall,  the 
pressure  rise  extended  from  0.14  in.  to  0.22  in.  The  first  pres- 
sure rise  indicated  by  the  center  trav- 
erse spread  between  0.56  in.  and  0.72 
in.  In  the  next  angle  shock  at  the  wall 
where  the  boundary  layer  was  thicker, 
the  pressure  rise  extended  from  0.80 
in.  to  0.95  in.  This  would  indicate  that 
angle  shocks  emanating  from  a  channel 
wall  are  unstably  anchored  to  a  sub- 
sonic boundary  layer  the  same  as  plane 
shocks. 

Shocks  at  High  Temperatures 

When  a  gas  flows  through  an  orifice  or 
a  parallel-walled  nozzle  under  a  pressure 
ratio  gieater  than  that  needed  to  pro- 
duce sonic  velocity,  the  sudden  adjust- 
ment of  the  gas  pressure  to  the  lower 
back  pressure  produces  a  variety  of  angle 
shock  patterns.  Fig.  8  shows  static  and 
impact  pressure  center  traverses  for 
such  flow  from  a  simple  nozzle,  to- 
gether with  the  probable  shock  and 
rare- faction  wave  pattern. 

The  fact  that  each  pressure  rise  is 
spread  over  a  distance  of  0.15  in.  indi- 
cates that,  owing  to  the  boundary  layer 
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in  the  nozzle,  the  sliock  pattern  was  not  stable  and  stationary 
in  space.  Further  evidence  of  such  instability  was  furnished  by 
the  whistling  of  the  impact  tube  in  the  region  from  0.05  in. 
inside  the  nozzle  to  0.10  in.  outside  of  it. 

The  loss  of  impact  pressure  and  the  subsequent  recovery  of  it 
between  0.15  and  0.40  in.  marks  the  region  where  the  flow  was 
sufficiently  supersonic  to  produce  a  serious  shock  in  front  of 
the  impact  tube.  The  distance  of  this  shock  bow  in  front  of  the 
impact  tube  appears  to  have  been  about  0.08  in. 

It  is  usually  assumed  that  the  first  angle  compression  shock 
emanates  from  the  nozzle  wall,  but  the  traverse  of  the  jet  from 
an  orifice,  shown  in  Fig.  9,  indicates  that  angle  compression 
shocks  can  originate  from  a  jet  boundary.  Sonic  velocity  was 
not  reached  imtil  0.2  in.  from  the  orifice  and  no  disturbance 
reached  the  center  traversing  tube  in  the  first  0.33  in.  At  that 
point  an  angle  shock  and  the  subsequent  rarefaction  were  fol- 
lowed by  a  plane  or  bow  shock  in  front  of  the  tube  shown  in  the 
stream. 

Such  shock  patterns  as  shown  by  the  curves  in  Figs.  8  and  9, 
respectively,  are  quite  common,  and  alternate  compressions  and 
rarefactions  usually  persist  for  6  or  more  oscillations  before  they 
are  damped  out.  With  low-temperature  gases,  they  are  detecta- 
ble only  by  making  pressure  traverses  or  by  schlerin  photo- 
graphs. 

When  these  shock  patterns  oc(;ur  in  high-temperature  prod- 
ucts of  combustion,  the  high-pressure-shock  diamonds  are  directly 
visible  and  are  usually  a  beautiful  blue  color.  Such  a  view  is 
illustrated  in  Fig.  10.  Fig.  11  shows  test  results  for  thermo- 
couple readings  taken  in  such  a  shock  region  in  a  stream  of 


jiroducts  from  an  orifice  approximately  2  in.  diam.  The  first 
reaction  to  such  results  as  shown  in  Fig.  11  is  that  of  utter  dis- 
belief. However,  the  chance  of  unusual  instrumentation  errors 
has  been  systematically  eliminated  by  repeating  such  test  runs 
under  a  variety  of  conditions,  always  with  the  same  results. 

The  next  idea  is  that  there  is  unburned  fuel  being  discharged 
from  the  orifice,  which  burns  on  the  surface  of  the  thermocouple 
and  produces  a  local  temperature  rise.  The  chance  of  this  being 
the  case  is  very  small  for  two  reasons:  (a)  The  original  burner 
temperature  checks  a  heat  balance  on  the  air  and  fuel  supplied 
within  2  per  cent;  (6)  the  appearance  and  color  of  the  shock 
diamonds  are  not  at  all  altered  by  the  insertion  of  a  thermocouple. 

Another  possible  explanation  is  that  some  unstable  oxides  of 
nitrogen  are  formed  in  the  region  of  shock  where  there  are  rates 
of  deceleration  greater  than  2,000,000  g.  This  again  has  been 
eliminated  by  burning  propane  with  pure  oxygen  and  observing 
exactly  the  same  phenomenon. 

Fig.  12  indicates  that  it  is  a  phenomenon  definitely  associated 
with  shocks.  The  data  for  the  curves  in  Figs.  12  and  13,  re- 
spectively, were  taken  in  a  '/le-in.  jet  with  a  fairly  large  thermo- 
couple ( Vs-in.  sheath),  so  the  cooling  effect  was  quite  pronounced. 
As  a  result,  the  excess  temperatures  shown  are  much  smaller 
than  those  shown  for  the  large  jet  in  Fig.  10,  but  they  are  more 
analytical. 

Fig.  12  indicates  further  that  this  phenomenon  is  definitely 
geared  to  the  strength  of  the  compression  shock,  for  it  com- 
pletely disappears  at  orifice  or  burner  pressure  ratios  where 
shocks  are  not  present.  It  would  naturally  be  expected  that  an 
unshielded  thermocouple  in  a  small  stream  of  hot  gases  would 
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Fig.  8    Sudden  Expansion  Fbom  a  Simplb  Nozzli 


Fig.  9     Sudden  Expansion  From  an  Orifice 
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'iG.  10     High-Pkessure-Shock  Diamonds  in  High-Temperature 
Products  of  Combustion 
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Fig.  11 


Fig.  12  (Above)     Thermocouple    Readings   in   Shock   Diamonds 

Using  Combustion  Products  of  Propane  and  Air:     Effect  of 

Burner  Pressure  Ratio  on  Excess  Temperature  Readings 

Fig.  13  (Below)     Thermocouple  Readings   in   Shock   Diamonds 

Using  Combustion  Producers  of  Propane  and  Air:     Effect  of 

Burner    Temperature     on     Excess    Temperature     Readings, 

Burner  Pressure  Ratio  =  3.62 

read  from  100  to  200  deg  F  low  owing  to  radiation  losses  and 
velocity  errors,  and  that  is  what  happened  below  pressure  ratios 
of  1.9. 

Fig.  13  indicates  that  within  the  limits  of  accuracy  of  the  test 
data,  the  amount  of  excess  temperature  is  fairly  independent  of 
burner  temperature  above  1100  F.  No  test  data  are  available 
for  initial  temperatures  between  1100  F  and  80  F.  However, 
repeated  attempts  to  find  excess  thermocouple  readings  in 
shocks  at  a  temperature  at  80  F  failed. 

Since  the  thermocouple  metal  reaches  a  temperature  greater 

than  the  total  temperature  of 
the  gas  before  it  is  in  energy 
equilibrium,  it  is  evident  that 
it  is  being  bombarded  by  some- 

thing  besides  gas  molecules. 

The  violent  deceleration  of  the 
hot  gases  in  the  shock  could 
produce  such  a  high  concen- 
tration of  ions  that  their  bom- 
bardment and  neutralization 
at  the  thermocouple  metal 
surface  could  cause  the  metal 
to  reach  a  very  high  tempera- 
ture before  energy  equilibrium 
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is   reached.     Tliis   phenomenon  deserves  a  thorough   investiga- 
tion.' 

Combustion  and  Flow 

Combustion  of  air-fuel  mixtures  is  customarily  done  with  no 
appreciable  flow  velocity  and  too  little  is  known  of  the  mechanism 
involved  when  combustion  occurs  in  a  high-velocity  air  stream. 
The  most  commonly  discussed  case  of  combined  combustion  and 
flow  is  that  of  burning  in  a  constant-area  passage  which  is  such  a 
short  distance  that  wall  friction  forces  may  be  neglected. 

If  it  is  also  assumed  that  the  area  is  constant  and  the  flow  is 
steady,  Equation  [7]  gives  the  condition  for  constant  momentum 
per  second  as 

[37] 


Pi     +    Plfl^     =     P2     +    P2V2^. 

solving  for  the  increase  in  velocity  pressure 


PiVj' 

2 


PlVl' 

2 


Pi 


[38] 


This  says  that  the  gain  in  velocity  pressure  is  only  one  half  the 
drop  in  static  pressure.  It  was  to  find  out,  if  possible,  how  the 
gas  stream  goes  about  losing  the  rest  of  this  static  pressure  with- 
out wall  friction  that  the  following  tests  were  made: 

The  I'/g-in-diam  thin-walled  tube,  used  in  Figs.  14  and  15, 
was  lined  up  axially  in  the  air  stream  from  a  3-in-diam  nozzle. 
Hydrogen  for  heating  was  introduced  through  the  '/s-in-diam 
tube  shown.     By  regulating  the  size  and  arrangement  of  hy- 

3  Such  a  program  is  being  carried  on  by  the  Mechanical  Engi- 
neering and  Physics  Departments  of  the  Rensselaer  Polytechnic 
Institute. 


drogen  jets  and  the  hydrogen  pressure,  a  uniform  temperatun 
could  be  reached  across  the  entire  cross  section.  Static  anc 
total  pressures  was  obtained  with  Vio-in.  traversing  tubes,  anc 
temperatures  with  an  unshielded  thermocouple  of  the  sam( 
maximum  diameter.  All  values  shown  represent  averages  o 
the  annulus  around  the  hydrogen  tube. 

In  general,  two  quite  different  sets  of  conditions  could  be  ob 
tained.  Wlien  the  heating  was  uniform  across  the  tube  the  re 
suits  in  Fig.  14  are  representative  of  several  tests.  The  burning 
was  a  series  of  explosions  at  the  resonant  frequency  of  the  tube 
This  frequency  was  high  enough  for  manometers  and  thermo- 
couples to  give  steady  readings,  but  it  was  a  violently  noisy  forn 
of  combustion.     The  thing   which   characterizes   these   results 


showa  in  Fig.  14  is  that  the  velocity  pressure 


(f) 


mcreases  ir 


the  region  where  the  temperature  is  rising. 

However,  when  the  hydrogen  tube  was  placed  slightly  off  centci 
and  the  tube  wall  became  a  bit  hotter  on  one  side,  the  combustior 
became  anchored  and  was  very  quiet  and  orderly.  The  traverses 
in  Fig.  15  are  typical  of  these  tests,  and  they  are  characterized 
by  the  fact  that  practically  all  of  the  velocity-pressure  increase 
occurred  ahead  of  the  temperature  rise.  The  chief  loss  in  total 
pressure  occurred  in  the  burning  region  where  the  velocity  pres- 
sure remained  essentially  constant. 

This  would  indicate  some  form  of  jet  separation  in  the  tube 
ahead  of  the  combustion,  followed  by  combustion  at  essentially 
constant  velocity  pressure  as  the  flow  again  filled  the  tube, 
Fig.  16  shows  the  velocity  plotted  against  absolute  temperature 
and  for  the  tube  to  flow  full,  the  velocity  would  be  essentially 
proportional  to  temperature,  that  is,  if  the  weight  of  hydrogen 


CuRyt_Ho.i-3_: 


[ u  IE 


t  t 


1  a  i  1  s 


L. 


Fig.  14     Combustion  in  Straight  Tube;     Explosive  Burning 


Fig.   15     Combustion  in  Straight  Tube;      Quiet  Burnin(; 
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idded  and  the  small  pressure  change  are  neglected.  The  points 
)efore  and  after  combustion  do  lie  very  close  to  such  a  line,  but 
he  velocity  in  the  burning  region  is  much  higher.  This  furnishes 
I  good  indication  that  the  tube  was  not  flowing  full. 

In  the  second  section  of  the  paper  it  was  concluded  that  even 
vith  Unsteady  flow,  Equation  [33]  would  be  true  for  averaged 
eadings.  For  the  case  of  the  curve  in  Fig.  14,  for  intermittent 
iow  and  explosive  burning,  the  ratio  of  velocity-pressure  increase 
..o  static-pressure  decrease  was  0.65  instead  of  the  theoretical  0.50. 
fn  the  steady-flow  case  in  Fig.  15  it  was  0.57,  and  for  another 
series  of  eight  tests  it  varied  from  0.52  to  0.58.  Since  Equation 
|[38]  does  not  account  for  the  addition  of  hydrogen  to  the  flowing 
i;as  stream  this  is  not  a  bad  check. 

It  is  easy  to  assume  that  just  because  combustion  occurs  in  a 
ponstant-area  tube  the  burning  occurs  at  constant  area.  The 
imited  amount  of  test  work  at  hand  would  seem  to  cast  doubt 
jon  the  validity  of  the  simultaneous  assumptions  of  steady  flow 
[and  constant-area  combustion. 

Some  of  the  implications  of  these  assumptions  may  be  had  by 
assuming  a  rate  of  heat  release  of  H  Btu  per  sec  per  cu  ft  of  vol- 
ume, and  stating  that  for  steady  operation  an  amount  of  energy 
AHdx  must  be  transported  out  of  the  volume  Adx,  each  second; 
or,  from  the  energy  equation 


W 
HAdx  =  WCJT  +  ~  vdv. 

Jg 


[39] 


For  acceleration,  with  no  wall  friction 

dv 
AdP  =  — pAdx  —  =  — pAvdv  — 


-pAdx  —  = 
dt 


W 


For  steady  flow 


W  =  pAgv  =  -^^v. 


dv [40] 


141] 


Integrating  Equation  [40] 


W 


AP  =  —  —  V  +  Ci. 
9 


For  any  initial  condition  Pi  and  Vi,  Ci  may  be  evaluated  to  give 

P) [43] 


W 


From  Equations  [40]  and  [43] 

vdv  = ~.dP 

W 

Differentiating  Equation  [41] 


^^  +  t^P^-P) 


A  A  A 

dT  =  -—-diPv)  =  -—LvdP  +  —-Pdv. 

RW  RW  RW 


[44] 


[45] 


From  Equations  [40],  [43],  and  [45] 


'^'RW"' 


"^> 


1-P)    J- 


gA^PdP 
RW 


From  Equations  [44]  and  [46]  in  [39] . 

dP HW 

dx 


.^il|' 


Wvi 
lAg 


+  (Pi  —  P) 


Jip 


R 


Since,  from  Equation  [41] 


PiA 

Wvi  =  WiVi  = Vi' 

RTi 


[46] 


[47] 


[48] 


/- 

lA' 

/ 

Vt 

L-  iTtiAeE'^Ff.aTl.iSt  Im 


i->m::im>..\,::\f^'l:  \  m-^t-.i    i:,-:a'r::h  .:r^-f  t 


Fig.  16     Velocity-Temperature  Traverse  for  Curve,  Fig.  15 


[42]       and  from  definition 


ygRTi 


=  Ml' 


and 


C„ 


1       C„ 


1 


R        J        R       /(t  — 1) 
Equations  [48],  [49],  and  [50]  with  Equation  [47]  given 


dP 
dx 


HJW{y  —  1) 


PigA 


1 


(y  +  1)  --(^Mi^-hl) 
^1  _l 


[49] 


[50] 


[51] 


This  type  of  a  pressure  gradient  has  all  the  earmarks  of  in- 
stability, for  {dP/dx)  becomes  a  greater  negative  value  as  P  be- 
comes smaller.  This  would  indicate  that  the  effect  would  tend 
to  cumulate  until  (dP/dx)  becomes  minus  infinity  which  would 
be  a  discontinuous  front.  Consequently,  analysis  appears  to 
bear  out  the  experimental  conclusion  that  true  constant-area 
burning  cannot  also  be  a  steady-flow  phenomenon. 
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